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i ii Abstract: We update the study of the production of direct J/ip in association with a Z 

12 boson at the Next-to-Leading Order (NLO) in a s by evaluating both the yield differential 

13 in Pt and the J/tp polarisation in the QCD-based Colour-Singlet Model (CSM). Contrary 
u to an earlier claim, QCD corrections at small and mid Pt are small if one assumes that the 

^-j- 15 factorisation and the renormalisation scales are commensurate with the Z boson mass. As it 

16 can be anticipated, the t-channel gluon-exchange (t— CGE) topologies start to be dominant 

17 only for Pt >mz/2. The polarisation pattern is not altered by the QCD corrections. This 
^vq is is thus far the first quarkonium-production process where this is observed in the CSM. 

w Along the same lines, our predictions for direct T + Z are also given. 

> 
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Introduction 





36 A few years ago, non-perturbative effects associated with colour-octet (CO) channels [1— 

37 3] were considered to be the only plausible explanation for the numerous puzzles in the 

38 predictions of quarkonium-production rates at hadron colliders. The situation has slightly 

39 changed since then, with the first evaluations of the QCD corrections [4-8] to the yields 

40 of J/V> and T (commonly denoted Q hereafter) produced in high-energy hadron collisions 

41 via Colour-Singlet (CS) transitions [9]. It is now indeed widely accepted [10-12] that a* 

42 and corrections to the CSM are significantly larger than contributions at mid and 

43 large Pt and that they should be taken into account in any analysis of their Pt spectrum. 

44 Nowadays, it not clear anymore that CO channels dominate and they are the only source 

45 of quarkonia. As a result, there is no consensus on which mechanisms are effectively at 

46 work in quarkonium hadroproduction at high energies, that is at RHIC, at the Tevatron 

47 and, recently, at the LHC. 

48 Polarisation predictions for the CS channel are also strongly affected by QCD correc- 

49 tions as demonstrated in [6, 8, 13, 14]. At NLO, Q produced inclusively or in association 

50 with a photon are expected to be longitudinally polarised when Pt gets larger, whereas 

51 they were thought to be transversely polarised as predicted at LO in the CSM [15, 16]. 

52 Such a drastic change is understood by the dominance of new production topologies. This 
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53 also explains the significant enhancement in the production rates as observed for increasing 

54 Pt- 

55 The situation is rather different at low Pt, where the CS predictions for Q at LO [9] 

56 and NLO [4-6] accuracy are of the same magnitude at RHIC energies; this shows a good 

57 convergence of the perturbative series. They are also in agreement [17-19] with the existing 

58 data from RHIC [20] energy all the way up to that of the LHC [21-28]. CO channels are 

59 most likely not needed to account for low Pt data -and thus for the Pt integrated yields. 

60 This is at odds with earlier works, e.g. [29], which wrongly assumed that \c feed-down 

61 could be the dominant CSM contribution. This is supported further by the results of recent 

62 works [30] focusing on production at e + e~ colliders which have posed stringent constraints 

63 on the size of C = +1 CO contributions which can be involved in hadroproduction at low 

64 Pt- Finally, this is reminiscent of the broad fixed-target measurement survey of total cross 

65 sections [31] which challenged the universality of the CO MEs. 




(a) Born (b) NLO loop (c) NLO real emis- (d) i-CGE (e) t-CGE 

sion from the heavy 
quark 



Figure 1: Representative diagrams contributing to J/ip and T (denoted Q) hadropro- 
duction with a Z boson in the CSM by gluon fusion at orders aa 2 s (a), aa^ (b,c,d) and 
initiated by a light-quark gluon fusion at order aaf (e) . The quark and antiquark attached 
to the ellipsis are taken as on-shell and their relative velocity v is set to zero. 

66 In this paper, we focus on the production of J ftp (and T) in association with a Z boson. 

67 Whereas this process may give us complementary information on quarkonium production if 

68 it happens to be experimentally accessible at the LHC, it also offers an interesting theoret- 

69 ical playground for the understanding of the QCD corrections in quarkonium-production 

70 processes. Our motivation was twofold: first, to see if the polarisation pattern of the J /if} 

71 is altered by the QCD corrections at large Pt', second, to see how large the effect of new 

72 topologies opening at NLO is, by comparing a full NLO computation to a simplified one - 

73 NLO*- with a infrared (IR) cut-off and neglecting loops. Our attention has also been drawn 

74 to this process by a previous analysis of the yield at NLO [32] which showed an intriguing 

75 result where NLO corrections were large at low Pt and getting smaller at large(r) Pt- 

76 Such a result could only be explained by a negligible contribution from new kinematically 

77 enhanced topologies and a large (positive) contribution from loop corrections at low Pp. 

78 As we shall demonstrate, the conclusion drawn in [32] are misguided by an unconventional 

79 choice of the factorisation and renormalisation scales (fip and /j,r), -way below mz— and 

80 a Pt range not large enough -compared to mz- to be able to observe the dominance of 
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81 t— CGE topologies. As a matter of fact, if one chooses a value for the scales commensurate 

82 with mz, rather than the transverse mass of the J/tp as done in [32], the NLO corrections 

83 are found to be small at small Pt- On the other hand, for Pt > mz/2, the NLO corrections 

84 are enhanced by a kinematical factor P^. 

85 The paper is organised as follows. In sections 2 and 3, we describe the evaluation of 

86 the cross section at LO and NLO accuracy in the CSM. We also explain how the partial 

87 NLO* yield is evaluated. In section 4, we present our results which we first compare to 

88 those from [32] with the same scale choice, at the same energy and in the same kinematical 

89 region. Then, we show our predictions in an extended Pt range for hf,r commensurate 

90 with mz and we discuss the ratio NLO over LO. We also study the sensitivity of our 

91 prediction on the aforementioned scales. Afterward, we compare the NLO* yield with the 

92 full NLO and we comment on the dependence on the IR cut-off at large Pt and on the 

93 impact of the t— CGE topologies. In section 5, we analyse the yield polarisation at LO, 

94 NLO and NLO*. In section 6, we give and discuss our predictions for T. Section 7 gathers 

95 our conclusions. 



97 



2 Cross section at LO accuracy 

In the CSM [9], the matrix element to create a 3 Si quarkonium Q with a momentum Pq 

98 and a polarisation A accompanied by other partons, noted j, and a Z boson of momentum 

99 Pz is the product of the amplitude to create the corresponding heavy-quark pair, Ai(ab — > 

00 QQ), a spin projector A r (A|si,S2) and 22(0), the radial wave function at the origin in the 

01 configuration space, obtained from the leptonic width, namely 



M(ab^Q\P Q ) + Z(p z )+j)= 



N(X\ Sl ,s 2 ) S u ' 22(0) 



V'»Q VN c V4n (2 .i) 
x M(ab -> Q?Q$(p = 0) + Z( Pz )+j), 

102 where Pq = pq+Pq, p = (pq —pq)/2, si and s 2 are the heavy-quark spins, and 5 n> / '\fN~ c is 

103 the projector onto a CS state. 2V(A|si, s 2 ) can be written as s 2 )^^u(^ , s\) in 

104 the non-relativistic limit with being the polarisation vector of the quarkonium. Summing 

105 over the quark spin yields to traces which can be evaluated in a standard way. 

106 At LO, there is only a single partonic process at work, namely gg — > J/ipZ -completely 

107 analogous to gg — > J/iJyy for J/ -0-prompt photon associated production- with 4 Feynman 

108 graphs to be evaluated. One of them is drawn on Fig. 1 (a). The differential partonic cross 

109 section is readily obtained from the amplitude squared 1 , 

no from which one obtains the double differential cross section in Pt (Pt = Pj/^t) an d the 

in J/i/j rapidity, y, for pp — > J/ipZ after convolution with the gluon PDFs and a change of 



1 The momenta of the initial gluons, fci,2, are, as usual in the parton model, related to those of the colliding 
hadrons (pi,2) through ki,2 = £i,2Pi,2- One then defines the Mandelstam variables for the partonic system: 
s = sxix 2 , i = (ki - Pj/4,) 2 and u = (fc 2 - Pj/^) 2 - 
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variable: 



da 



dydP 7 



v /s( v /s-m T e-») : 



dx 



itit 



2sP T g(x 1} fi F )g(x 2 {x 1 ), fj, F ) da 



dt 



(2.3) 



114 3 Cross section at NLO accuracy 

us The NLO contributions can be divided in two sets: one gathers the virtual corrections 

116 which arise from loop diagrams, the other gathers the real (emission) corrections where 

117 one more particle appears in the final state. In the next sections, we briefly describe how 
us these are computed. 

119 3.1 Virtual corrections 

120 The computation of the virtual corrections involves three types of singularities: the ul- 

121 traviolet (UV), the infrared (IR) and the Coulomb ones. UV divergences arising from 

122 self-energy and triangle diagrams are cancelled after renormalisation. A similar renormal- 

123 isation scheme as in Ref. [33] is used, except for the fact that, in the present study, the 

124 bottom quark is also included in the renormalisation of the gluon field. The renormali- 

125 sation constants Z m , Zi and Z% which are associated to the charm quark mass m c , the 

126 charm-field if) c and the gluon field are defined in the on-mass-shell (OS) scheme while 

127 Z g , for the QCD gauge coupling constant a s is defined in the modified minimal-subtraction 

128 (MS) scheme: 
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6Z 
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(3.1) 



129 where is Euler's constant, (3q = — \^F n f is the one-loop coefficient of the QCD 

130 beta function and nt is the number of active quark flavours. We take the three light quarks 

131 u,d,s as massless and consider the quarks c and b as heavy; therefore nj=5. In SU(3) C , 

132 we have the following colour factor: Tp = \,Cp = \-,Ca = 3. Finally, f3' = (3$ + |Tp = 

133 -g-CU — ^Tpnif where nif = rif — 2 = 3 is the number of light quark flavours. 

134 After having fixed our renormalisation scheme, there are 111 virtual-correction dia- 

135 grams, including counter-term diagrams. Diagrams that have a virtual gluon line connect- 

136 ing the charm quark pair forming the J/ip lead to Coulomb singularity ~ 7r 2 /|p|, which 

137 can be isolated and mapped into the cc wave function. 
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138 The loop integration has been carried out thanks to the newly upgraded Feynman 

139 Diagram Calculation (FDC) package [34], with the implementation of the reduction method 
wo for loop integrals proposed in Ref. [35]. 

ui 3.2 Real corrections 

142 The real corrections arise from three parton level subprocesses: 

g + g -> J/^ + Z + g, (3.2) 
g + q(q) -> J/^ + Z + q(q), (3.3) 
q + q -> J/il> + Z + g, (3.4) 

143 where q denotes light quarks with different flavours (u,d,s). We have not considered the 

144 contributions from the processes cc^-J/ip + Z + g and g + c(c) — > J/ip + Z + c(c). Both 

145 are IR finite and can be safely separated out from the other ones. The charm- gluon fusion 

146 contribution may be non-negligible in the presence of intrinsic charm. It will be considered 

147 in a separate work. 

us The contribution from the quark-antiquark fusion (Eq. (3.4)) is also IR finite and small. 

149 The phase-space integration of the other two subprocesses will generate IR singularities, 

150 which are either soft or collinear and which can be conveniently isolated by slicing the 

151 phase space into different regions. We use the two-cutoff phase-space-slicing method [36], 

152 which introduces two small cutoffs to decompose the phase space into three parts. The 

153 real cross section can then be written as 

„Real „Soft , ,Hard Collinear , ,Hard NoncoUinear /o r\ 

a = a + a + a . 1,3. 5J 

154 The hard noncollinear part u Hard Noncollmear [ s tr finite and can be numerically com- 

155 puted using standard Monte-Carlo integration techniques. Only the real subprocess of 

156 Eq. (3.2) contains soft singularities. Collinear singularities appear in both real subpro- 

157 cesses of Eq. (3.2) and Eq. (3.3), but only as initial-state collinear singularities. As shown 

158 in Ref. [36], all these singularities can be factored out analytically in the corresponding 

159 regions. When combined with the IR singularities appearing in the virtual corrections (see 

160 section 3.1), the soft singularities of the real part cancel. Yet, some collinear singularities 

161 remain. These are fully absorbed into the redefinition of the parton distribution function 

162 (PDF): this is usually referred to as the mass factorisation [37]. All the singularities are 

163 thus eventually analytically cancelled. 

164 3.3 NLO* cross section 

165 In order to evaluate the NLO* contributions, we use the framework described in [38] based 

166 on the tree-level matrix-element generator MADONIA [39] slightly tuned to implement an 

167 IR cut-off on all light parton-pair invariant mass. The LO cross section has also been 
lea checked with MADONIA. 

169 The procedure used here to evaluate the leading- Pt NLO contributions is exactly the 

170 same as in [8] but for the process pp — > J ftp + Z + jet. Namely, the real-emission contri- 

171 butions at aa^ are evaluated using MADONIA by imposing a lower bound on the invariant 
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172 mass of any light parton pair (s™ m ). The underlying idea in the inclusive 2 case was that 

173 for the new channels opening up at NLO which have a leading- Pt behaviour (for instance 

174 the i-CGE), the cut-off dependence should decrease for increasing Pt since no collinear or 

175 soft divergences can appear there. For other NLO channels, whose Born contribution is 

176 at LO, the cut would produce logarithms of •Sij/s™ 111 , which are not necessarily negligible. 

177 Nevertheless, they can be factorised over their corresponding Born contribution, which 

178 scales as P^ S , and are thence suppressed by at least two powers of Pt with respect to the 

179 leading- Pt contributions (P^ 6 ) at this order. The sensitivity on s™ m should vanish at large 

180 Pt- This argument has been checked in the inclusive case for T [8] and ip [10] as well as 

181 in association with a photon [14]. Because of the presence of the Z boson mass, it is not a 

182 priori obvious that t— CGE topologies dominate over the LO ones. It is thus not clear at all 

183 how such procedure to evaluate the NLO* yield can provide a reliable evaluation of the full 

184 NLO of J/ip + Z. In fact, at mid Pt, significantly below the Z boson mass, the difference 

185 of the Pt dependence of the NLO and LO cross sections is maybe not large enough for the 
lee dependence on s™ m to decrease fast. Having at hand a full NLO computation, we can carry 

187 out such a comparison and better investigate the effect of QCD corrections in quarkonium 

188 production. This is done after our complete results are presented. 

189 4 Results for J/tp + Z: differential cross section in Pt 

190 4.1 Comparison with Mao et al. [32] 

191 In order to compare our results with those of [32], we take Vs = 14 TeV and \y J ^\ < 

192 3.0. We also set the factorisation and renormalisation scales at the same value, namely 

193 fjt F = hr = m J T ^ = yjm*^ + pj. We also take a = 1/137, |i?j/^(0)| 2 = 0.91 GeV 3 , 

194 m c = 1.5 GeV, m z = 91.1876 GeV and sin 2 (#vi/) = 0.23116. Our LO results (also cross- 

195 checked with MADONIA) do match with those of [32] (compare both blue curves on Fig. 2). 

196 However, as depicted in Fig. 2, we are not able to reproduce the NLO results presented 

197 in the later reference. At low Pt, we have found a K factor smaller than one {i.e. the yield 

198 at NLO is smaller than at LO) while they obtained a value larger than one. The way a s 

199 is precisely evaluated in both computations may differ but this can hardly explain a sign 

200 change in the aa^ contributions. There is also a difference in the renormalisation scheme: 

201 we have included both charm and bottom quarks in the renormalisation of the gluon field 

202 contrary to what has been done in the previous analysis. Yet, we do not believe that this 

203 could explain the discrepancies between both results. 

204 That being said, a scale close to mz, rather than the transverse mass of the J/ip 

205 taken in [32], seems more appropriate as done for instance for Z + b— jet [40]. This has an 

206 important effect on the scale sensitivity, less on the final numbers predicted for the yields, 

207 as we shall discuss in the next section. 

208 In Fig. 3, we show the scale sensitivity at low Pt around two different choices of the 

209 "default" scale value, /Jo, (a) the transverse mass of the J/ip and (b) the Z boson mass. We 

2 "Inclusive" is used here in opposition to "in association with another detected particle" which is indeed 
a more exclusive process. 
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Figure 2: Comparison between our results (solid lines) and that of Mao et al. [32] (dashed 
lines) for the differential cross section for J/tp + Z vs. the J /ip Pt at LO (blue) and NLO 
(gray) with [ip = nn = m 
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Figure 3: (a) Renormalisation and factorisation scale dependence of the LO and NLO 



yield for Pt > 3 GeV with /iq 



m 



J/i> 



. (b) Same plot as (a) for fiQ = mz- 



210 emphasise that we believe the latter choice to be more appropriate owing to the presence 

211 of the Z boson in the hard process. One sees that around mz (b), the cross section at 

212 NLO is more stable, except for the bump at 2mz which can be corrected by properly 

213 setting the value of in the running of coupling constant (currently 0.151 MeV with 
2u mt = 180 GeV), which matters for [Ir > mt- The NLO results are clearly unstable at low 

215 scales and they may then artificially be enhanced. In the following sections, we investigate 

216 further the dependence of the scale sensitivity for different domains of the J/ip transverse 

217 momenta. 
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218 4.2 Results for the differential cross section in Pt at y/s = 8 TeV and 14 TeV 

219 In the following, we show our results for \y J ^\ < 2.4 -the usual J/ip acceptance for the CMS 

220 and ATLAS detectors- at 8 TeV and 3 14 TeV and for the renormalisation and factorisation 

221 scales set at mz- We have kept 4 the cut P^ > 3 GeV. 

222 The parameters entering the cross-section evaluation have been taken as follows: \Rj/^(0)\ 2 

223 0.91 GeV 3 , Br(J/> l+£~) = 0.0594, m c = 1.5 GeV with m J/v , = 2m c , m b = 4.75 GeV, 

224 a = 1/128. Our result at ^/s = 14 TeV are depicted on Fig. 4. The dotted blue line is our 

225 LO result and the solid gray line is our prediction at NLO. It is obvious, contrary to what 

226 was obtained in [32], that the yield at NLO is getting larger than at LO for increasing 

227 Pt- This is similar to what happens in the inclusive case. This is also indicative that new 

228 leading Pt topologies, in particular t-CGE, start dominating rather early in Pt despite of 

229 the presence of a Z boson in the process. At Pt = 150 GeV, the NLO yield is already ten 

230 times that of LO. 
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Figure 4: Differential cross section for J ftp + Z vs. Pt at \fs = 14 TeV at LO (blue 
dashed) and NLO (gray solid) with = f^R = n^z along with the NLO* for different 
values of s™ m (red dotted, yellow double dotted, purple dash dotted and pink long-dash 
dotted). 



231 The dominance of t-CGE topologies can be quantified by a comparison with the results 

232 from the NLO* evaluation. As aforementioned, because of the Z boson mass, it was not a 

233 priori clear that the NLO* evaluation could make any sense here. Indeed, as long as the 

234 contribution from the sub-leading Pt topologies are significant, the NLO* would strongly 

3 The cross section at 13 TeV is 12 % smaller than at 14 TeV. 

4 Note that we could have evaluated the cross section for lower Pt where the cross section is well behaved. 
However, we do not expect -at least in the central region- any experimental measurement to be carried 
out in this region owing to the momentum cut on the muons because of the strong magnetic fields in the 
ATLAS and CMS detectors. 
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235 depend on the arbitrary IR cutoff 5 which is used to mimic the effect of the loop contributions 

236 which regulate the soft gluon emission divergences. We are in a position to check from which 

237 Pt the NLO* starts to reproduce the full NLO and becomes to be less sensitive on the IR 

238 CUt. 

239 The various dotted lines on Fig. 4 show the NLO* evaluation for different cut-off values. 

240 Two observations can be made: 1) they converge to the NLO steadily for increasing Pt, 2) 

241 for Pt > mz, the NLO* evaluations are within a factor of 2 compatible with the complete 

242 NLO yield. This confirms that loop corrections are sub- leading in Pt and can be safely 

243 neglected for Pt larger than all the masses relevant for the process under consideration 

244 and that new topologies appearing at NLO, the i-CGE ones, dominate at large Pt- At low 

245 Pt, where the NLO and LO yield are similar, the NLO* overestimate the NLO. 




Figure 5: Differential cross section for J/ip + Z vs. Pt at y/s = 8 TeV at LO (blue dashed) 
and NLO (gray solid) with up = = mz- 

246 As regards the possibility to study such a process at the LHC, the Pt differential cross 

247 sections times the branching Br(J/?/; — > at the smallest Pt accessible by ATLAS and 

248 CMS (3 to 5 GeV depending on the rapidity) is of the order of 1 fb/GeV at 14 TeV (Fig. 4) 

249 and three times less at 8 TeV (Fig. 5). These do not take into account the branching of the 

250 Z in dimuons (~ 3%), but these are central values which can be affected by a factor of 2-3 

251 of theoretical uncertainties. In the most optimistic case, by integrating on the accessible 

252 Pt range, by using both muon and electron decay channels for the J/ip, by expecting an 

253 indirect cross section of 40 % and by detecting the Z boson with hadronic channels such 

254 that it could be detected 40 % of the time, it may be envisioned to detect something like 

255 four hundred events at 14 TeV with 100 fb _1 of data -possibly a factor 2-3 above given 

256 the theory uncertainty. At 8 TeV with the 20 fb _1 of data expected to be collected in 

257 2012, we expect only about thirty events to be recorded. Clearly, there are more promising 

258 processes, such as J/tp + 7 [13, 14] or J/ip + D [17, 41], to learn more on the production 

5 Not to be confused with the cutoff used in the full NLO computation, on which the final results does 
not depend. 
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259 mechanisms of the J ftp- Nevertheless, the study of J ftp + Z may suffer less from trigger 

260 limitations and could thus still be at reach at the LHC. In any case, it is an ideal theory 

261 playground to analyse the effects of QCD corrections on quarkonium production, which 

262 have been the key subject in the recent years in the field. 

263 4.3 Scale sensitivity at different Pp 

264 From the observations made above, we expect the real emission contributions at aa^ to 

265 dominate for Pt > mz/2. This should therefore impact on the scale dependence of the 

266 yield. At low Pt (<C mz), we expect a reduced scale dependence since we really deal with a 

267 process at NLO accuracy. At large Pt, the leading process is pp — > J/ip + Z + parton. The 

268 loop contributions are not expected to reduce the scale sensitivity since they are small. 

269 On the contrary, we expect a larger sensitivity on the renormalisation scale, hr, since 

270 the leading process shows an additional power of q^/Ur). In practice, we study the scale 

271 sensitivity by varying \ip and \ir together and then [ir alone by a factor 2 about the 

272 "default" scale mz with 3 cuts in Pt -i.e. 3, 50 and 150 GeV. 




(a) Scale dependence at LO 



1 

u/u, 

(b) Scale dependence at NLO 



Figure 6: Scale dependence of the yield at LO (a) and NLO (b) for P T > 3 GeV, Pt > 
50 GeV, Pt > 150 GeV where both the renormalisation and factorisation scales are varied 
together (fip = fiR, solid lines) about fio = mz and only the renormalisation scale is varied 
(fip fixed, dashed lines). Note that a has been kept fixed. 



273 On Fig. 6, we do observe, as anticipated for PT^rnz (red curves), a stronger scale 

274 sensitivity of the NLO yield (b) -at aa^- than of the LO yields (a)- at cto? s . The NLO 

275 curve with \xp fixed clearly shows that the sensitiviy essentially comes from fiR. At mid Pt 

276 (orange curves), the scale sensitivities are similar at LO and NLO, while at low Pt (black 

277 curves), the NLO yield is less scale dependent than the LO -in agreement with the common 

278 wisdom regarding the NLO computations. Note also that the 3 LO curves showing the sole 

279 dependence of /ir are identical since one can factor out a common a 2 s since our choices of 

280 fiR do not depend on Pt- 
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281 5 Polarisation: polar anisotropy in the helicity frame 



282 The polar anisotropy of the dilepton decay of the J/ip, \e or a, can be evaluated from the 

283 polarised hadronic cross sections: 

284 To evaluate a(Pr), the polarisation of J/ip must of course be kept throughout the 

285 calculation. The partonic differential cross section for a polarised J/ip is expressed as: 

^ = a e(A) • e*(A) + £ a„ • e(A) Pj • e*(A), (5.2) 
ij=l,2 

286 where A = Ti,T2,L. e(Ti), e(T2), e(L) are respectively the two transverse and the longi- 

287 tudinal polarisation vectors of J ftp', the polarisations of all the other particles are summed 

288 over in n dimensions. One can find that a and aij are finite when the virtual and real 

289 corrections are properly handled as aforementioned. There is therefore no difference in 

290 the partonic differential cross section d&x/dt whether the polarisation of J/ip is summed 

291 over in 4 or n dimensions. Thus, we can just treat the polarisation vectors of J/ip in 4 

292 dimensions. There are usually several different choices of the polarisation frames, as dis- 

293 cussed in Ref. [42-44]. In our calculation, we have chosen to work in the helicity frame. 

294 The polarisation can be obtained in a given frame by taking the corresponding polarisation 

295 vectors in Eq. (5.2). 
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Figure 7: (a) Pt dependence of the polarisation (or azimuthal anisotropy) in the helicity 
frame of the direct J/tp produced with a Z boson at LO, NLO and NLO* (for 2 values of 
the IR cut-off) at ^/s = 14 TeV. (b) Same as (a) at LO and NLO at = 8 TeV. 



296 Our results at 14 TeV in Fig. 7 (a) clearly show that the direct- J/tp yield in association 

297 with a Z boson is increasingly longitudinally polarised in the helicity frame for increasing 



P. 



J/i> 



The NLO and NLO* results coincide and the latter is nearly insensitive to the 
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299 IR cutoff. Interestingly, the NLO and the LO results are also very similar. This is the 

300 first time that such a robustness of the polarisation against QCD corrections is observed 

301 for the colour-singlet channels. For the J/ijj produced inclusively or in association with a 

302 photon, the yield at LO and NLO are found to have a completely different polarisation. 

303 Our interpretation is that, when a Z boson is emitted by one of the charm quarks forming 

304 the J ftp, the latter is longitudinally polarised, irrespective of the off-shellness and of the 

305 transverse momentum of the gluons producing the charm-quark pair. This is not so when 

306 a photon or a gluon is emitted in the final state. In the present case, we also note that the 

307 polarisation at 8 TeV (Fig. 7 (b)) is nearly exactly the same as at 14 TeV. 

308 6 Results for T + Z 

309 Along the same lines as for J ftp, we have also evaluated the cross section and the polar- 

310 isation for direct-T production in association with a Z boson. Experimentally, the CDF 

311 Collaboration at Fermilab has set a 95 % C.L. upper value for such a cross section at 

312 y/8 = 1.8 TeV [45], namely 

a(pp -»• T + Z + X) x Br(T -> < 2.5 pb. (6.1) 

313 Further studies with the entire data set recorded by CDF is under process [46]. At 
3u y/s = 1.8 TeV, a quick evaluation of the total cross section (without y cut, nor Pt cut) 

315 gives, for the CSM, a value close 6 to 0.1 fb (~ 0.2 fb by taking into account a similar 

316 feed-down fraction (~ 50%) than that for the inclusive case). A similar evaluation for 

317 the CO transitions is highly dependent on the choosen LDME values and on the expected 

318 impact of the feed-down. Values span from ~ 0.06 fb for the direct yield with CO LDMEs 

319 fit [47] from the early prompt Tevatron data, up to ~ 3.75 fb as evaluated in [48], passing 

320 by ~ 0.4 fb for the direct yield using CO LDMEs fit from the latest Tevatron results taking 

321 into account some NLO QCD corrections [49]. 
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Figure 8: Differential cross section for direct T + Z vs. Pt at LO (blue-dashed) and NLO 
(gray-solid) with fi F = fi R = m z at 14 TeV (a), 8 TeV (b) and 7 1.96 TeV (c). 



We should however keep in mind that these are central values and that theoretical uncertainties can be 
of the order of a factor 2-3. 
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322 This, is in any case, significantly below the CDF upper bound obtained with 83 pb 

323 of data. Given these small theoretical values, we fear that such process cannot be experi- 

324 mentally accessed at the Tevatron, unless contributions from colour-octet transitions, from 

325 double-parton interactions or from feed-downs are unexpectedly large. At the LHC at 14 

326 TeV, the expected yield in the CSM for the central rapidity region accessible by CMS and 

327 ATLAS is of the order 5 fb (still including the branching of the T in muons). The central 

328 values for the differential cross sections vs. Pt at LO and NLO are shown in Fig. 8 (a-c). 

329 An enhancement by a factor 2 to 4 can certainly be expected if the feed-downs from excited 

330 bottomonium states and the usual theoretical uncertainties are taken into account. 

331 By comparing Fig. 8 (a-c), one also notices an interesting phenomenon: the NLO and 

332 LO yields start to depart from each other at low Pt for increasing y/s. This can probably 

333 be attributed to an increasing -negative- size of the loop corrections in this region at small 

334 x. This is in fact reminiscent to what has been observed in the inclusive case [4, 6, 17, 18]. 

335 In the latter case, the situation is worse since the NLO cross section can become negative 

336 for large yfs and small Pt- It hints at significant NNLO corrections at small Pt in the 

337 small- x region, as anticipated in [50]. 

338 For the sake of the comparison with the J/tp case, we have also computed the polari- 

339 sation at LO and NLO. As it can be seen on Fig. 9, the yield polarisation at LO and NLO 

340 are very alike, though slightly different from for the J/tp case -most probably due to the 

341 change in the quarkonium mass compared to the Z mass. 

1 
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Figure 9: Pt dependence of the polarisation (or azimuthal anisotropy) in the helicity 
frame of the direct T produced with a Z boson at LO and NLO at y/s = 14 TeV. 

7 We have considered a wider rapidity range than usual for the CDF quarkonium analyses since CMX 
muons can be used in such a correlation analysis owing to the smaller background compared to inclusive 
measurements [46]. 
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342 7 Conclusions 



343 In conclusion, we have studied the effects of the QCD corrections to the production of direct 

344 J/ip and T via colour-singlet transitions in association with a Z boson at the LHC. We 

345 have found, contrary to an earlier study [32], that the NLO QCD corrections are consistent 

346 with the expectations, namely increasing for increasing Pt and small at low Pt- We expect 

347 that a few hundred J ftp + Z events could be detected at the LHC at 14 TeV with 100 fb _1 

348 of data. At 8 TeV with 20 fb — 1 , there may be just enough events to derive a cross section 

349 and not only an upper bound on its value. Interestingly, the CSM yield expected for direct 

350 T + Z is of the same order of magnitude than that of direct J ftp + Z at 14 TeV, if not 

351 larger. 

352 We have studied the scale sensitivity of the J/ip + Z cross section at LO and NLO. At 

353 low Pt, it is smaller when QCD corrections are taken into account. On the contrary, at large 

354 Pt, i.e. when Pt > mx/2, the dominant contributions are of the kind of gg — > Q + Z + jet. 

355 These involve an additional power of a s and the sensitivity on the renormalisation scale is 

356 larger. Along the same lines of arguments, one expects a further increase of the CSM cross 

357 section at large Pt when the leading- Pt contribution contained in the NNLO topologies 

358 are accounted for. That being said, the presence of the Z boson mass renders the CSM 

359 prediction more precise at low Pt compared to the inclusive case, for which the leading- Pt 

360 contributions at NLO and NNLO dominate at lower Pt- 

361 We have also found that the yield polarisation is not altered by the QCD corrections. 

362 From this observation, we have concluded that when a Z boson is emitted by one of the 

363 heavy quarks forming the quarkonium, both the J/ip and the T are longitudinally polarised 

364 at LO and NLO, thus independently of the off-shellness and of the transverse momentum 

365 of the gluons producing the heavy-quark pair. This is at odds with the cases of inclusive Q 

366 production and Q + 7 production, and this motivates further theoretical and experimental 

367 investigations. 
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